On the spin asymmetry of ground states in trapped two-component Fermi gases with 

repulsive interactions 
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We examine the spin asymmetry of ground states for two-dimensional, harmonically trapped 
two-component gases of fermionic atoms at zero temperature with weakly repulsive short range 
interactions. Our main result is that, in contrast to the three-dimensional case, in two dimensions 
a non-trivial spin-asymmetric phase can only be caused by shell structure. A simple, qualitative 
description is given in terms of an approximate single particle model, comparing well to the standard 
results of Hartree-Fock or direct diagonalization methods. 

PACS numbers: 05.30Fk, 03.65.Ge 



I. INTRODUCTION 

Cold atomic gases with multiple components have been 
realized for example in Bose-Einstein condensates 
01 • The nowadays possible experimental studies of atomic 
Fermi gases [5, 6, 7, 8, a, E El, [H, [ll have raised 
much recent interest in the different underlying quantum 
phases, with a main attention on pairing. Fermi gases 
with spin imbalance have been a hot topic more recently 
in connection with superfluidity issues While most 



of the previous studies focus on attractive interactions, 
Duine and MacDonald pointed out that the repul- 
sive regime might offer the possibility of ferromagnetism 
in ultra-cold atom gases with, as they say, 'unprecedented 
experimental control'. Describing spin (magnetic) prop- 
erties of Fermi systems is an important issue in many- 
body physics. The trapped dilute fermionic atom gases, 
with their flexibility in tuning important system param- 
eters such as the interaction strength, could provide op- 
portunities to look for new magnetic phenomena and to 
test the theoretical models of magnetic properties. 

Possible spin symmetry breaking for a Fermi gas in 
the unitarity limit [3, E3, [3, [3, EO] was discussed by 



Chevy (2l|. Maruyama and Bertsch [22| have included 



time-dependent, external magnetic flelds, yielding spin- 
asymmetric ground states. 

Sogo and Yabu jisi studied the three-dimensional, 
trapped Fermi gas with two components. (The more 
general multi-component problem was discussed earlier 
by Salasnich et al. [2J|)- These components can be 



characterized by different intrinsic spin states, or for ex- 
ample, two hyperflne states, playing the role of pseu- 
dospin in a two-component Fermi gas. For weak inter- 
actions between the components, the zero-temperature 
ground state is symmetric, with equal up- and down- 
densities and vanishing total spin or pseudospin. Sogo 
and Yabu showed, that in three spatial dimensions, de- 
pending on the interaction strength between the compo- 
nents, a spin-asymmetric, so-called "collective ferromag- 
netic" state may become energetically favorable. This 
spin or pseudospin asymmetry in the ground state arises 



as it becomes energetically favorable to increase the 
Fermi energy of one component rather than increasing 
the interaction energy between the two species. Sogo and 
Yabu [iS] based their investigation on the Thomas-Fermi 
approximation 
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26l |. an eflicient way to compute the 
ground state properties of a large system - however, at 
the cost of ignoring the shell corrections [ii*]. 

These studies were performed for a three-dimensional 
harmonic trap. However, by appropriate manipulation 
of the trapping potential, conflned ultra-cold atoms can 
achieve conditions of reduced dimensionality. We thus 
found it useful to further investigate the question of spin 
asymmetry for an isotropic harmonic trap in two dimen- 
sions. (In one dimension, a uniform fermi system can be 
solved exactly with a Bethe ansatz [27|. Recently, this 
was also generalized to a harmonic trap within the local 
density approximation [j^). 

We found that in two dimensions the spin asymmetry 
can only arise from the shell-structure and Hund's rule. 



II. THE THOMAS-FERMI APPROXIMATION 
FOR A TWO-COMPONENT FERMI GAS 

Consider fermionic particles of equal mass m at zero 
temperature, with two internal (spin) states labeled by 
cr e {u,d}. The fermions shall be subject to repulsive 
two-body contact interactions. In atomic units (h — Lo — 
m = 1), our model Hamiltonian in two dimensions is 
(r=|r|) 
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where Ng- denotes the number of particles in each of 
the two spin states. Here, only the repulsive, weak- 
interaction regime of Eq. ^ is considered. The cou- 
pling constant g oc a/lz > is proportional to the 
ratio of the three-dimensional scattering length a and 
the length scale Iz of the ground state in the frozen-out 
z— dimension. From the condition that the (mean-fleld) 
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interaction energy should be much smaller than the Fermi 
energy 

gnl|^lF « 1, (2) 

with n° denoting the particle density (of one of the 
species) in the center of the trap, one is led to the di- 
luteness condition 



< g < 27r. 



(3) 



We start by treating the problem within the Thomas 
Fermi (TF) approximation which should be valid in the 
limit of large particle numbers where the specific contri- 
bution from shell structure plays a minor role. In the 
Thomas-Fermi approximation, the total energy is 

1 



E 



dr 



(4) 



where n„ and and n = n„ + are the (pseudo)spin- 
and total densities, respectively. The Euler-Lagrange 
equations for the densities, while fulfilling the constraint 
of constant particle number N — + Nd = J dr n, are 
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In order for the system to be in equilibrium we must have 
Mtt = Mrf = M- Further on for the system |[5]) to be a stable 
solution we must have 
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> ^ g < 27r, 



(6) 



which is guaranteed by the diluteness condition. By view- 
ing (lU as a matrix equation in the variables n„ and 
it has the system determinant 47r^ — g^- So for g < 27r, 
which is again guaranteed by the diluteness condition, 
the system ^ has a unique algebraic solution which cor- 
responds to a stable minimizer of Eq. (|4]) 



nd 



n,, 



(7) 



This TF result suggests that there should be no sta- 
ble spin-asymmetric density configurations in the two- 
dimensional system under study. This is qualitatively dif- 
ferent from the corresponding TF investigation for three- 
dimensional system where a non-trivial phase separation 
called "collective ferromagnetic state" has been reported 
[23, 2^- This is due to the different exponent of the 
kinetic energy functional in the three-dimensional case, 
leading to non-linear equations for the densities. 

Inserting the density n = Uu + rid = 2n„, into the 
energy functional the total energy equals 
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(8) 



We then note that for g = 2i: this energy coincides with 
the fully polarized phase, i.e. 1/2 times Eq. ^ with 
g = and N — > 27V, hence giving a physical motivation 
for the condition in Eq. ([6]). This describes the trivial 
phase separation to a totally polarized system similar to 
what can happen in three dimensions [23j. 



III. THOMAS-FERMI SINGLE-PARTICLE 
MODEL 



Let us now incorporate the shell structure. Obviously, 
the density profile given by the above Eq. ^ will also be 
valid for two unequally populated non-interacting {g = 0) 
spin states, as according to Eq. ([5]), n„ and are then 
uncoupled. This is the situation when particles of one 
(pseudo)spin state are filling a higher shell than the other 
spin state, such that \^iu ~ Md| — ^eff {^eff is defined in 
Eq. ifTT]) ). The major effect of the repulsive interaction 
is to redistribute atoms from the center of the cloud to 
larger r values. This is qualitatively modeled by the pro- 
file of the upper Eq. |[7| also when — Mdl — ^eff 
(and agrees well with a numerical solution of Eq. ^ for 
g <^ 27r). Moreover, the profile of Eq. ([7]) does not break 
the SU{2) symmetry, as also was shown to be the case 
for the exact density (29| . 

We now begin from the single-particle Hartree-Fock 
equations [i^l for the Hamiltonian ^ in the weakly re- 
pulsive regime 



-^A + gnd (x, y) + i + 



Ei^uf%,u- (9) 



Now, for the purpose of presenting a simple analytical 
model, we approximate the particle density with the g- 
dependent TF density corresponding to Eq. ([7]) 

Ud ix,y) « (kl) = 



1 



2TT + g 



where fid = \/(2 H- g/Tr) Nd and {■) is Heaviside's step- 
function. After separating the wavefunctions ipi^u {x,y) 
in Eq. ^ into the Cartesian coordinates we obtain e.g. 
for the x-dimension 



1l'^,u {x) = (10) 




g^^dO {fid - x'^/2) 



2tt- 
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Figure 1: (Color online) The Nu-Nd chart for harmonically trapped 2D fermions with two spin states according to TFSPM 
(here with g — 0.1). The black dots (•) show the states with lowest energy for a given N — Nu + Nd, the corresponding spin 
for this states are seen in the inset figure. The green triangles (A) show case ii) {spin symmetry) from equation l|16p . green 
circles (o) show case in) (polarized). The red stars (★) are the magic numbers Mj = j (j + 1) = 2, 6, 12, .... 



By re-arranging the terms in Eq. (fTTj) , this can be viewed 
as a one-dimensional harmonic oscillator with the pertur- 
bation 

We then apply first order non-degenerate perturbation 
theory in each Cartesian dimension 

£i,u ^ i + + {4>i\hp\(j)t) , (13) 

where is the eigenfunction i — 0, 1, 2, ... to the one- 
dimensional harmonic oscillator with frequency uj = 1. If 
we now only keep the leading order terms in 5 <C 27r, we 
can for practical purposes write Eq. (fT3l) as 

e^,u~^+l + ^ f (^2N'd-^xA\cj),ix)\''dx. 

(14) 

Clearly this is the x-dependent part of what one obtains 
if one treats the mean-field interaction term in Eq. ([9]) 



as gnuo from the beginning, where nno refers to the 
TF-density of the non-interacting system. 

Obviously, Eq. (fT3|l takes the value for the non- 
interacting HO (-i.u = « + 1/2 in the two limits g ^ 
and TVrf — + (Nu —> N). In total, we have two Cartesian 
coordinates and two spin-components such that the total 
energy is |25|] 

Etot ^ Eu+Edfii2'^ (ei^u + ei,d) - g dr n„nd. (15) 

occ. •' 

At zero temperature we have a Fermi-Dirac occupation of 
each component such that the HO main quantum num- 
ber increases by one unit whenever a shell is filled. There 
are three special cases for which the total energy (fTSl can 
easily be written down explicitly. 

i) The non-interacting case: When g = we have 
e^,^ = i + l/2. 

ii) Spin symmetry: When Nu — Nd, before subtract- 
ing the last term in Eq. ifTSj) . the total energy is sim- 
ply twice the one for a single component system with 
the number of particles being N/2 and with a poten- 
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tial V{g,N/2) = uj'lgr'^ [2 + Cshtft, but in this case 
the Thomas-Fermi radius Rtf — V^Ptf is equal for 
both components. Then all the single particle wavefunc- 
tions decay very rapidly for r > Rtf and we can for- 
mally drop the factors 9 (fid — 1"^) in Eq. (flTI) . then 

Ei.rr = We^gP (i + 1/2) + tshift- 

Hi) The totally polarized case: When N = and 
Nd = (or vice versa) the interaction disappears which 
gives Ci^u = i + 1/2. The cases i), ii) and Hi) are all 
captured by the formula js^l 



EuM {g, a) 



27r ri , /Vl + 8a- r 

-ceil I X 

27r + 5[6 V 2 



6a + 1 



ceil 




where ceil (•) is the closest integer from above. For case i) 
one has Etot = (5 = 0, a = Nu) + Ed{g ^ 0,a ^ Nd). 
In case ii) Etot = 2i?„^d (5, " = -^/2) — e where e — 
g fdr n2/4 = gN^/"^ / ^ySGn'^ + ISng. In case Hi) Etot = 
Eu,d ig^O,a = N). 

For all other cases, Eq. (fT4|) can be calculated (numer- 
ically) for each single particle level Ci^a which are then 
summed up according to Eq. lfT5|) . If one compares this 
procedure with a case where Eq. (flGj) is applicable e.g. 
case ii) where Nd — Nu — N/2, one finds that the rel- 
ative difference between those two approaches is small 
(e.g. < 2 • 10"'^ for g — 0.1 and N = 10) and decreases 
with TV. 

The approximative model we suggest here to illus- 
trate shell-induced spin asymmetry relies upon the fact 
that a spherically symmetric TF particle-density utf oc 
(/I — r^/2) does not break the SU (2) symmetry, which 
is present for the non-interacting problem due to the 
harmonic trap. It has been shown in lowest order per- 
turbation theory that the ^-degeneracy for a given HO 
shell is not lifted, even when the exact expression for the 
particle-density is used [1^. This fact supports the va- 
lidity of our approximative "Thomas-Fermi single particle 
model", TFSPM, for weak interactions. However, except 
for the noninteracting case {g = 0), the TFSPM obvi- 
ously overestimates the many-body energy compared to 
more sophisticated computational methods, where there 
is more freedom to vary the many-body wavefunction. 
TFSPM successfully finds the correct spin-configurations 
with lowest energy, see Figure [U and reproduces the shell 
correction with very high accuracy, see Figure [2l What 
we here refer to as the correct spin-configuration is what 
we obtained from full Hartree-Fock calculations (see be- 
low) for all presented values of N. For < 12 we con- 
firmed the same spin-configuration by direct numerical 
diagonalization of the Hamiltonian |(T|) . 



IV. COMPARISON TO THE HARTREE-FOCK 
METHOD 

Now we turn our interest to the full Hartree-Fock 
calculation. We calculate the particle densities (e.g.) 
nd{x,y) = J2i'=i\fi,df (without demanding spherical 
symmetry) for each spin component iteratively from the 
mean-field single-particle states (pi^d [x,y) corresponding 
to the Hartree-Fock Eq. 

For each total number of particles = Nu+Nd, we de- 
termine the lowest-energy state self-consistently. For this 
state, the total spin is defined as S* = (iV„ — Nd) /2. A 
zig-zag pattern with increasing amplitude and frequency 
occurs, see inset of Figure [J, which is identical for the 
Hartree-Fock calculations and the TFSPM. This pattern 
(inset of Figure 1) is a direct consequence of Hund's rule. 
The trapped atoms at a degenerate shell can lower their 
interaction energy by maximizing the number of atoms 
of the same species - in other words, by aligning their 
spins. This mechanism in contact-interacting fermionic 
systems leads to Hund's first rule and what we here call 
shell-induced spin asymmetry, in close similarity to long- 
range interacting systems. Here, however, Hund's rule 
has a more dramatic effect since it removes completely 
the interaction between atoms of the same spin state. 

However, in absolute values, the energy difference as- 
sociated with this pattern is a decreasing function of the 
number of particles. To be more specific, the lowest en- 
ergy state carrying the largest total spin in each shell is 
Smax (j) ={j + l)/2, where j = 1, 2, 3, ... is the index 
of the shell. Since the magic numbers of an isotropic 
two-dimensional harmonic oscillator with two spin-states 
are M, = 2^^-^^(771+1) = j (j + 1) - 2, 6, 12, ... , 
Smax is related to the total number of particles as 
Smax — VN /2. The energy difference of the state with 
S {N) = Smax compared to the state S (N) = is of the 
order AE ~ gVN. Comparing this to the leading order 
term in the TF estimate of the total energy of the state 
with S{N) = from Eq. © E 2N^/^/3, the relative 
energy gain of the shell-induced spin asymmetry scales 
as AE/E ~ g/N, which becomes irrelevant for large TV. 
This is in agreement with the analysis in section II, which 
had shown that the Thomas-Fermi approximation could 
not predict any stable spin-asymmetric phase for the sys- 
tem, but only the trivial, totally polarised system. 



V. CONCLUSIONS 

In summary, the short-range weakly repulsive inter- 
acting fermi gas is a simple toy model that illustrates 
very directly Hund's rule: The total spin is maximized 
at mid-shell, in order to minimize the total energy, what 
we here have refered to as shell-induced spin asymmetry. 
For the case of an isotropic two-dimensional harmonic 
trap we have illustrated a simple approximative single 
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Figure 2: (Color online) Oscillating part of the total energy 
(shell correction) for the state with lowest energy, at a given TV 
(here with <; = 0.1). The solid black line with stars (*) shows 
the TFSPM, the red dashed line with circles (o) the numerical 
Hartree-Fock results. The green triangles (A) show case ii) 
{spin symmetry) from equation l|16p. which is only defined for 
even A''. In this case the total energy was ~ 0.3% lower for 
Hartree-Fock compared to TFSPM and the deviations in the 
shell correction are very small. Investigating even smaller val- 
ues of g the curves approach each other further. The smooth 
part of the energy is constructed such that the different curves 
for the oscillating part should coincide at the minima of the 
shell correction (i.e. at the magic numbers) and they follow 
a line y oc A''^'^^. 

particle model (TFSPM), that qualitatively reproduces 
the energy landscape of the system, when compared to 
Hartree-Fock calculations. This model is restricted to 
two dimensions and weak interactions since it relies upon 
the fact that the TF density profile does not break the 
SU (2) symmetry of the non-interacting problem and the 
fact that perturbation theory is used. This also outrules 
the possibility to see super-shell structure, reported in 
the three-dimensional case [3l|, within this approach. 

ACKNOWLEDGEMENTS 

This work was financially supported by the Swedish 
Research Council and the Swedish Foundation for Strate- 
gic Research, as well as the European Community project 
ULTRA-ID (NMP4-CT-2003-505457). 



[1] M.H. Anderson, J.R. Ensher, M. R. Matthews, C.E. Wie- 
man and E.A. Cornell, Science 269, 189 (1995). 



[2] K.B. Davis, M.O. Mewes, M.R. Andrews, N.J. van 
Druten, D.S. Drufee, D.M. Kurn and W. Ketterle, Phys. 
Rev. Lett. 75, 3969 (1995). 

[3] W. Ketterle, Rev. Mod. Phys. 74, 1131 (2002). 

[4] E. A. Cornell and C. E. Wieman, Rev. Mod. Phys. 74, 
875 (2002). 

[5] B. De Marco and D.S. Jin, Science 285, 1703 (1999). 
[6] M.J. Holland, B. DeMarco, and D.S. Jin, Phys. Rev. A 

61, 053610 (2000). 
[7] S.R. Granade, N.E. Gehm, K.M. O'Hara, and J.E. 

Thomas, Phys. Rev. Lett. 88, 120405 (2002). 
[8] S. Jochim, M. Bartenstein, G. Hendl, J. Hecker Den- 

schlag, and R. Grimm, Phys. Rev. Lett. 89, 273202 

(2002) . 

[9] C.A. Regal and D.S. Jin, Phys. Rev. Lett. 90, 230404 

(2003) . 

[10] Z. Hadzibabic, S. Gupta, C. A. Stan, C. H. Schunck, 

M. W. Zwierlein, K. Dieckmann, and W. Ketterle, Phys. 

Rev. Lett. 91, 160401 (2003). 
[11] M. Greiner, C.A. Regal and D.S. Jin, Nature (London) 

426, 537 (2003). 
[12] M. W. Zwierlein, C. A. Stan, C. H. Schunck, S. M. F. 

Raupach, S. Gupta, Z. Hadzibabic, and W. Ketterle, 

Phys. Rev. Lett. 91, 250401 (2003). 
[13] C.A. Regal, M. Greiner and D.S. Jin, Phys. Rev. Lett. 

92, 040403 (2004). 
[14] M. W. Zwierlein, A. Schirotzek, C.H. Schunck, and W. 

Ketterle, Science 311, 492 (2006). 
[15] R.A. Duine and A.H. MacDonald, Phys. Rev. Lett. 95, 

230403 (2005). 

[16] "The Many-Body Challenge Problem'^ formulated by G. 

F. Bertsch, 1999, see R. F. Bishop, Int. J. Mod. Phys. B 

15, iii (2001). 
[17] K. M. O'Hara et al.. Science 298, 2179 (2002); 
[18] J. Kinast, A. Turlapov, J. E. Thomas, Q. Chen, J. Stajic, 

and K. Levin, Science 307, 1296 (2005); 
[19] T.-L. Ho, Phys. Rev. Lett. 92, 090402 (2004); 
[20] A. Bulgac, J.E. Drut, and P. Magierski, Phys. Rev. Lett. 

96, 090404 (2006). 
[21] F. Chevy, " Universal phase diagram of a strongly in- 
teracting Fermi gas with unbalanced spin populations", 

cond-mat/0605751 (2006). 
[22] T. Marujama and G.F. Bertsch, Phys. Rev. A 73, 013610 

(2006). 

[23] T. Sogo and H. Yabu, Phys. Rev. A 66, 043611 (2002). 

[24] L. Salasnich, B. Pozzi, A. Parola and L. Reatto, J. Phys. 
B: At. Mol. Opt. Phys. 33, 3943 (2000). 

[25] M. Brack and R. K. Bhaduri: Semiclassical Physics 
(Westview Press, Boulder, USA, 2003). 

[26] P. Ring and P. Schuck, The Nuclear Many-Body Prob- 
lem, (Springer- Verlag, New York, 1980). 

[27] M. Gaudin, Phys. Lett. 24A, 55 (1967). 

[28] H. Hui, L. Xia-Ji and P.D. Drummond, "FFLO 
phase in one-dimensional polarized Fermi gases", 
cond-mat/0610448 

[29] B. P. van Zyl, R.K. Bhaduri, A. Suzuki, and M. Brack, 
Phys. Rev. A 67, 023609 (2003). 

[30] M. Ogren, unpublished (2004). 

[31] Y. Yu, M. Ogren, S. Aberg, S.M. Reimann and M. Brack, 
Phys. Rev. A 72, 051602(R) (2005). 



